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This  project  deals  with  optimal  synthesis  of  an  RC  high-pass 
ladder  network.  It  shows  how  to  obtain  a  network  with  the  maximum 
possible  gain  and  minimum  total  series  capacitance.  The  current 
transfer  function,  source  and  load  conductances  are  specified.  Lagrange 
multipliers  method  is  used  to  obtain  optimization.  The  entire  process 
is  completed  using  a  microcomputer  which  can  handle  a  large  number 
of  poles  and  can  perform  the  optimization  in  a  very  short  time. 
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I.  INTRODUCTION 

Due  to  our  increasing  dependence  on  integrated  circuit 
technology,  optimal  synthesis  of  filters  has  become  more 
and  more  important.  The  size  of  the  chip  is  directly 
proportional  to  the  amount  of  resistance  and  capacitance 
contained  within.  Thus,  by  minimizing  these  values,  one 
can  minimize  the  cost  of  each  chip.  This  method,  however, 
does  not  necessarily  reduce  the  overall  cost  of  production. 

The  man-hours  involved  in  arriving  at  these  optimum  values 
may  more  than  offset  the  savings  realized  by  smaller  chips. 
This  is  especially  true  as  the  complexity  of  the  filter 
increases.  To  alleviate  this  problem,  one  can  use  a  micro¬ 
computer  to  aid  in  design.  This  reduction  in  man-hours 
makes  optimal  synthesis  economically  desirable. 

The  optimal  synthesis  of  a  double-terminated  high-pass 
RC  filter  is  very  tedious  if  done  by  hand.  The  whole  process, 
however,  can  be  done  in  only  a  fraction  of  the  time  using  a 
microcomputer  which  has  some  form  of  BASIC.  The  program  is 
limited  only  by  the  storage  capacity  and  accuracy  of  the 
microcomputer. 

II.  SYNTHESIS  PROCEDURES  AND  MICROCOMPUTER  APPLICATION 


To  synthesize  a  double-terminated  high-pass  RC  filter,  as 
depicted  in  Fig.  1,  one  must  first  specify  a  high-pass  current 
transfer  function  of  the  form: 
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where  the  poles  are  simple  and  on  the  negative  real  axis.  The 
source  and  load  conductances  must  also  be  specified.  From  the 
current  transfer  function  one  can  find  the  open-circuit  impedance. 
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This  then  can  be  expanded  into  partial  fractions 
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The  residues  of  each  pole  can  now  be  determined. 
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The  calculation  of  the  residues  of  z^2  takes  place  in  lines 
350-430  of  the  example  program  shown  in  the  Appendix,  and  are 
contained  in  the  list  A(I). 

If  the  network  is  compact  then: 
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where  k^  are  the  residues  of  z^  and  k^  are  the  residues 
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From  the  previously  defined  quantities  z^  has  the  form: 
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These  two  equations  must  satisfy  the  following  constraints. 


By  using  these  constraints  and  the  Lagrange  multipliers 
method,  one  can  obtain  the  following  quadratic  equation: 


Once  X  is  known,  it  can  be  used  to  obtain  the  residues  of  Zjj 
since  x^  *  Xou  .  The  value  of  X  and  the  residues  of  z^  are 
calculated  in  lines  440-490  of  the  example  program  with  C3 
containing  X  and  the  X(I)  containing  the  residues  of  z^. 

By  letting  s  =  0  and  substituting  the  expression  for  X 
into  one  can  obtain  an  equation  for  H. 


Once  the  residues  of  z^  and  H  have  been  determined,  one 
must  return  the  partial  fraction  expansion  of  to  a  single 
fraction  before  a  continued  fraction  expansion  can  be  performed. 
To  do  this  requires  that  one  be  able  to  expand  a  polynomial 
from  its  factored  form.  Consider  the  function. 
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where  the  coefficients  a.  must  be  determined.  It  can  be 
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shown  that  the  coefficients  a-  are  a  sum  of  products  following 
the  form: 

al  ‘  P1  *  p2  *  ***  *  Pn 

a2  “  plp2  +  plp3  *  +  pn-lpn 

a3  =  P1P2P3  +  P1P2P4  *  +  pn- 2pn- lpn 

an  "  Pi  '  ?2  ‘  «V"  Pn 

In  evaluating  these  coefficients  one  must  ensure  that  the 
number  of  FOR-NEXT  loops  is  independent  of  the  number  of 
poles  to  maintain  program  flexibility. 

To  do  this,  consider  the  following  truth  table: 

p3  p2  pl 

0  0  0 

0  0  1 

0  1  0 

0  1  1 

10  0 
1  0  1 

1  1  0 

1  1  1 

A  FOR-NEXT  loop  can  be  used  to  form  the  product  of  Pj  and  pn 
for  each  row.  When  the  pole  is  low,  however,  it  is  skipped 
and  the  counter  is  not  incremented.  The  product  is  then  summed 
into  a  list  where  the  position  is  determined  by  the  value 
of  the  counter  when  the  product  is  completed.  Thus  F(l) 
corresponds  to  a^,  F(2)  corresponds  to  and  so  forth.  To 
generate  the  truth  table  one  can  use  a  sinusoid  with  the 
frequency  depending  on  the  particular  column.  When  the 
sinusoid  is  negative,  the  pole  is  considered  high.  When  it 
is  positive  the  pole  is  low.  Lines  520-810  perform  the 
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function  of  combining  the  partial  fraction  expansion. 
Zjj  is  recombined,  one  can  use^n  algorithm  to  obtain  ; 
2  expansion  of  the  admittance  function.  The  algorithm 
shown  below. 
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The  element  values  are  given  by: 
a. 


a 


11 

21 
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Line  980-1020  perform  the  Cauer  2  continued  fraction  expansion 
using  the  array  S(J,I). 

Finally,  the  element  values  are  printed  in  a  convenient 
form  along  with  the  network  gain. 


III.  EXAMPLES 

A.  Consider  the  following  current  transfer  function  with  a 
load  conductance  of  5  mhos  and  a  source  conductance  of  1  mho. 

*o  Hs3 

17  17-T-rrTF-^TrT^T^ 

Solving  for  2 ^  ,  yields: 
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Knowing  the  residues  of  2^-,,  one  can  solve  for  X  and  H, 
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With  X  known  one  can  now  determine  the  residues  of  2^  and 


thus  the  expression  for  2^ 
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can  now  be  expanded  into  a  continued  fraction  using  an 
algorithm  for  Cauer  2  as  shown  in  Table  I. 

TABLE  I 


1 

2.351 

5.694 

.3148 

11.53 

.5734 

5 


6 

6 

2.552 

.4482 

1.424 

.1235 

.2153 

.0431 


11 

6 

1 

8.448 

2.  598 

.0431 

3.402 

.9569 

.  3475 

.0431 

.7118 

.0431 

The  completed  circuit  is  shown  in  Figure  2. 

The  entire  synthesis  process  can  be  done  using  the  micro¬ 
computer  program  shown  in  the  Appendix.  One  must  enter  3  for 
the  number  of  poles,  1  for  the  source  conductance,  5  for  the 
load  conductance,  and  each  of  the  poles  1,  2,  3  entered  sep¬ 
arately.  From  this  data  the  program  prints  out  the  value  for 
each  circuit  element  and  the  gain  of  the  circuit. 

B.  Consider  the  following  current  transfer  function  with  a 
load  and  source  conductance  of  1  mho. 
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Knowing  the  residues  of  z12»  one  can  solve  for  X  and  H 

1(2  *  1777^"  x  -  1  *  0 
X  =  1 


H  =  1 


1  +  1277 


=  7.303  x  10 


-2 


With  X  known  one  can  now  determine  the  resiudes  of  z ^  which 
in  this  case  are  the  absolute  value  of  the  residues  of  z^* 


*»  — 

07303  £ 

•  „  115.7  x  10'6  ^ 

.  2709  A  13 

.18  L  92.47  .  134.6" 

Z11 
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+  21.58s4  +  678. 

Is3  +  6537s2  +  20280s  +  14400 

Z11 

s5 

+  55s4  +  1023s3 

+  7645s2 

+  21076s  +  14400 

Yn  can 

now  be 

expanded  into  a 

continued 

fraction  using  an 

algorithm  for  1 

Cauer  2  as  shown 

in  Table 

II 

TABLE 

II 

14400 

21076  7645 

1023 

55  1 

1 

14400 

20281  6537 

678.1 

21.58  .0730 

18.11 

795.3 

1108  344.9 

33.42 

.9269 

3.525 

225.6 

292.6  72.93 

4.792 

.0730 

2.957 

76.32 

87.85  16.53 

.6695 

2.322 

32.86 

24.05  2.813 

.0730 

1.027 

23.00 

10.00  .5000 

2.322 

13.  78 

2.300  .0730 

2.957 

4.660 

.3304 

3.524 

1.322 

.0730 

18.11 

.0730 

1 

.0730 

The  completed  circuit  is  shown  in  Figure  3, 


Once  again,  the  process  can  he  done  using  a  microcomputer 
in  a  similar  fashion  for  Example  A. 

IV.  CONCLUSION 

As  evidenced  by  the  preceding  examples,  the  use  of 
microcomputers  greatly  reduces  the  time  required  to  optimally 
synthesise  a  filter.  One  is  required  to  know  only  the  source 
conductance,  load  conductance  and  each  pole  to  obtain  the 
circuit  gain  and  the  value  of  each  circuit  element.  In 
addition  to  decreased  time,  the  microcomputer  also  provides 
a  greater  degree  of  reliability  over  long-hand  methods. 

This  is  expecially  true  as  the  filter  becomes  more  complex. 

Reduced  time  and  greater  reliability  makes  the  use  of 
a  microcomputer  essential  in  the  optimal  design  of  filters. 


.. 
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100 
110 
120 
130 
140 
150 
lbO 
170 
180 
190 
200 
210 
220 
230 
240 
2  50 
260 
2  70 
2  80 
290 
300 
310 
320 
330 
340 
34  2 
344 
350 
360 
370 
380 
390 
400 
410 
420 
4  30 
440 
450 
460 
470 
480 
490 
492 
500 
502 
504 
510 
520 
530 


PRINT  "WOULD  YOU  L 1  Ki-  A  DESCRIPTION  OF  THE  PROCRAM?" 

PRINT  "ENTER  1  FOR  YES,  0  FOR  NO" 

INPUT  V 

IF  V-0  THEN  220 
DIM  R(ll)  ,F(11) ,T(21) 

DIM  S( 22,11) 

PRINT  "THIS  PROCRAM  SYNTHESIZES  AN  OPTIMAL 
PRINT  HIGH-PASS  FILTER.  FROM  THE  CURRENT 
PRINT  TRANSFER  FUNCTION  ONE  MUST  ENTER  EACH 
PRINT  POLE,  THE  SOURCE  CONDUCTANCE,  AND  THE 
PRINT  LOAD  CONDUCTANCE.  (THE  POLES  ARE 
PRINT  ENTERED  AS  POSITIVE  NUMBERS.)" 

LET  E-C-K-0 
LET  B-l 

LET  P-3. 141592645 

PRINT  "INPUT  THE  NUMBER  OF  POLES  (1<N<10)" 

INPUT  N 

PRINT  "INPUT  THE  SOURCE  CONDUCTANCE" 

INPUT  Cl 

PRINT  "INPUT  THE  LOAD  CONDUCTANCE" 

INPUT  G2 
FOR  l-l  TO  N 
PRINT  "INPUT  P" ; I 
INPUT  LCD 
NEXT  I 

REM  THE  FOLLOWING  STEPS  CALCULATE  THE  RESIDUCES  OF  Z12  WHICH 

REM  ARE  THEN  USED  TO  CALCULATE  THE  RESIDUES  OF  Zll 

FOR  1-1  TO  N 

FOR  .1-1  TO  N 

IF  .1-1  THEN  390 

LET  B-B*L(I)/*(L(J)-L(I)) 

NEXT  .1 

LET  A ( I ) -ABS  (B*L(I)) 

LET  B-l 

LET  C-OA(I)/L(I) 

NEXT  I 

LET  Cl - ( 1 -G2/G1 ) /C 

LET  C2-SQR(C1*C1*4*G2/G1) 

LET  C3-iC2-Cl)/2 
FOR  1-1  TO  N 
LET  X ( I ) -A( I ) *C3 
NEXT  I 

REM  WITH  THE  RESIDUES  OF  Zll  KNOWN,  ONE  CAN  CALCULATE  THE.  CAIN  H 
LET  H-G2/G1/ (1*C3*C) 

REM  THE  FOLLOWING  STEPS  ARE  USED  TO  PUT  THE  PARTIAL  FRACTION 

REM  EXPANSION  OF  Zll  BACK  INTO  ONE  FRACTION 

LET  D«2*N 

FOR  J«1  TO  D-l 

FOR  1-1  TO  N 
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540  IF  SIN(P*(J+.001)*2/2AI)>0  THEN  570 
550  LET  E  =  E  ♦  1 
560  LET  B=B*L ( I) 

570  NEXT  I 

580  LET  F(E)-F(E)+B 

590  LET  E  =  0  i 

600  LET  B-l 

610  NEXT  J 

620  LET  l)=l)/2 

630  FOR  J»1  TO  N 

640  FOR  1  =  1  TO  N - 1 

650  I<>J  THEN  670 

660  LET  K=K+ 1 

670  LET  K-K+l 

680  LET  M(I)=L(K) 

690  NEXT  I 
700  LET  K=0 
710  FOR  0=1  TO  D-l 
720  FOR  1=1  TO  N-l 

730  IF  SlN(P*(0+0.001)*2/2  I)  0  THEN  760 
740  LET  E  -  E  ♦  1 
750  LET  B-B*M(t) 

760  NEXT  I 

770  LET  Q(J,E)-Q(J,E)*B*X(J) 

780  LET  E-0 
790  LET  B*  l 
800  NEXT  0 
810  NEXT  J 
820  LET  R(l)-F(l) 

830  FOR  1=1  TO  N 
840  LET  R(1)-RU)*X(I) 

850  NEXT  I 

860  FOR  J-l  TO  N 

870  FOR  1=1  TO  N 

880  LET  R(J+1)-R(J+1)*Q(I,J) 

890  NEXT  I 

900  LET  R(J*1)-R(J+1)+F(J*1) 

910  NEXT  J 

920  LET  S(1,N*1)-1 

930  LET  S( 2  ,N+ 1 ) -H/G2 

940  FOR  1-1  TO  N 

950  LET  S(1,N*1-1)-F(I) 

960  LET  S(2,N*1-I)-R(I)*H/G2 
970  NEXT  I 

972  REM  FINALLY  A  CAUER  2  EXPANSION  IS  PONE 
980  FOR  J-l  TO  2*N 
990  FOR  1-1  TO  N 

1000  LET  S(J  +  2,I)  =  (S(J*l,l)*S(Jtl*l)-S(J,n*SlJM,l*n)/SlJM,l) 
1010  NEXT  I 
1020  NEXT  J 

1030  PRINT  "THE  ELEMENT  VALUES  ARE:" 

1040  FOR  1-1  TO  2*N*1 

1050  LET  T(t)-S(I,l)/S(I*l,l) 


1060  LF.T  U-2* INT (1/2) 

10  70  IF  II- 1  THEN  1110 

1080  LET  K-K+l 

1090  PRINT  "G" ; K ;T( I ) 

1100  GO  TO  1130 

1110  LET  E-E  +  l 

1120  PRINT  ”C";E;1/T(I) 

1130  NEXT  I 
1140  PRINT 

1150  PR  TNT  "THE  GAIN  OF  THIS  FILTER  IS";H 
1160  END 

Example  A 

FILTER 

WOULD  YOU  LIKE  A  DESCRIPTION  OF  THE  PROGRAM? 
ENTER  1  FOR  YES,  0  FOR  NO 
?  1 


THIS  PROGRAM  SYNTHESIZES  AN  OPTIMAL  HIGH- 
PASS  FILTER.  FROM  THE  CURRENT  TRANSFER  FUNCTION 
ONE  MUST  ENTER  EACH  POLE,  THE  SOURCE  CONDUCTANCE, 
AND  THE  LOAD  CONDUCTANCE.  (THE  POLES  ARE 
ENTERED  AS  POSITIVE  NUMBERS . 1 
?  3 

INPUT  THE  SOURCE  CONDUCTANCE 
?  1 

INPUT  THE  LOAD  CONDUCTANCE 
?  5 

INPUT  P  1 
?  1 

INPUT  P  2 
?  2 

INPUT  P  3 
?  3 

THE  ELEMENT  VALUES  ARE: 

G  1  -  1. 

C  1  -  0.425307 
G  2  -  5.69406 
C  2  -  3.17637 
G  3  -  11.5297 
C  3  -  1.74378 
C  4  -  5. 

THE  GAIN  OF  THIS  FILTER  IS  0.215297 


FILTER 


WOULD  YOU  LIKE  A  DESCRIPTION  OF  THE  PROGRAM 
ENTER  1  FOR  YES,  0  FOR  NO 
?  0 

INPUT  THE  NUMBER  OF  POLES  (1<N<10) 

?  5 

INPUT  THE  SOURCE  CONDUCTANCE 
?  1 

INPUT  THE  LOAD  CONDUCTANCE 
?  1 

INPUT  P  I 
?  1 

INPUT  P  2 
?  4 

INPUT  P  3 
?  9 

INPUT  P  4 
?  16 

INPUT  P  5 
?  25 

THE  ELEMENT  VALUES  ARE: 

G  1  -  1. 

C  1  =■  5. 52264E-2 
G  2  =■  3.  52454 
C  2  -  0.338232 
G  3  -  2.3223 
C  3  >  0.973783 
G  4  *  2.32244 
C  4  -  0.338229 
G  5  *  3.52454 
C  5  *  5 . 52264E- 2 
G  6  *  1. 

THE  GAIN  OF  THIS  FILTER  IS  7.30256E-2 


Figure  1 
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